Anomalous Behavior of the Spin Susceptibility of Strongly Correlated Fermi Systems 
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The spin susceptibility x{T) of strongly correlated Fermi systems is investigated in the density 
region where Fermi-liquid theory fails. We attribute this failure to a specific quantum phase tran- 
sition associated with a rearrangement of the Landau state at low temperatures T, retaining the 
assumption that the Landau quasiparticle picture survives in a generic sense. Taking into account 
the resulting modification of the quasiparticle distribution function, the spin susceptibility x(T) is 
shown to contain a Curie-Weiss component xcw(T) ~ (T — Ow)" 1 , with the Weiss temperature 
6w vanishing at the critical density for the transition. 
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The manifestation of non-Fermi-liquid behavior in 
strongly correlated Fermi systems provides valuable clues 
to a fundamental microscopic understanding of these sys- 
tems. This letter was stimulated by the findings of recent 
studies of laboratory simulants of two-dimensional (2D) 
liquid 3 He and the electron g aSl L2iMi£i£iL& Two dramatic 
effects have been revealed by data for the spin suscepti- 
bility x(T) and specific heat of 2D liquid 3 He and for the 
temperature dependence of Shubnikov-de Haas conduc- 
tivity oscillations in the 2D electron gas. First, the results 
point to a divergence of the effective mass M* in these 
systems at a critical density p^ where, concomitantly, 
the single-particle spectrum e(p) becomes flat. Second, 
as seen in the experimental data 1 ' 2,5 on %(T — ► 0), Lan- 
dau theory begins to fail even before this critical density 
is reached. 

Ordinarily, the non- Fermi- liquid behavior of \(T — > 0) 
is attributed to a disorder-driven localization transition. 
To be sure, disorder plays a part, as evidenced in dataS 
on the spin susceptibility of submonolayer 3 He at den- 
sity p — 0.0064A~ 2 (~ 0.1 layer), absorbed onto a Nuclc- 
pore substrate pre-plated with several 4 He layers. The 
spin susceptibility of this low-density disordered sys- 
tem is found to contain a significant Curie component 
Xc(T) r~j T , which progressively diminishes with in- 
crease of 4 He coverage D4, vanishing at D4 ~ 2.5 lay- 
ers. However, such a Curie component xc{T) is also 
found in the spin susceptibility of a high-mobility 2D 
layer in silicon^ In this case, disorder effects are sup- 
pressed, while the effective interaction between electrons 
is strong, since the r s parameter attains values around 
10. Furthermore, the same behavior of x(T — > 0) is 
observed 1 ^ when two monolayers of 3 He are absorbed 
onto a graphite substrate pre-plated with a monolayer 
of 4 He. The density of liquid 3 He in this experiment is 
about an order of magnitude higher than in the work of 
Ref.0. Again the effective interaction is strong, since M* 
is drastically enhanced. Thus, at least two examples ex- 
ist in which the interaction appears to be the dominant 
playmaker, rather than disorder. 

If indeed disorder proves to be irrelevant in some 



strongly interacting Fermi systems, how does the interac- 
tion give rise the inferred singular behavior of \(T — ► 0)? 
A possible answer to this question is provided by a 
microscopic mechanism proposed over a decade ago^S 
which predicted a flattening of the single-particle spec- 
trum in strongly correlated Fermi systems well before 
the phenomenon was observed in angle-resolved photoc- 
mission studies In this scenario, known as fermion 
condensationjiiiii^iSiiLiSiiS the flattening effect is as- 
sociated with a phase transition in which the conven- 
tional Landau state suffers a rearrangement. It is in- 
trinsic to this model that the Landau-Migdal quasiparti- 
cle picture retains its validity beyond the critical point. 
However, at T — the Fermi step np(p) — 9(pp — p) is 
replaced by a new momentum distribution, rio(p), deter- 
mined by the variational condition 



SE/5n(p) = n 



(1) 



The distribution no(p) differs from np{p) in a finite mo- 
mentum interval pi < p < pf that includes the Fermi mo- 
mentum pp. Since the variational derivative 5E/Sn(p) is 
by definition the energy e(p) of a Landau quasiparticle, 
the condition Q implies that in this domain, the single- 
particle spectrum must be completely flat, pinned to the 
chemical potential fi: 



e(p) = M > Pi<P<Pf ■ 



(2) 



The set of single-particle states having e(p) = p is called 
the fermion condensate (FC) in analogy with the Bose 
condensate existing in liquid 4 He, since both condensates 
entail similar (5-like singular contributions to the density 
of states. 

At finite T, the degeneracy of the FC spectrum is lifted 
in accordance with the formula 14 

^(p,T^0) = e(p)- A1 = Tlni^M , Pl <p< Pf , (3) 

which ensures consistency between the Fermi-Dirac ex- 
pression n(p,T) = [1 + exp(£(p)/T)] _1 and the solution 



n (p) of Eq. CJ at T = 0. At higher T the FC density 
falls to disappear at the critical temperature Tt. 

Having set the stage, let us now evaluate the spin sus- 
ceptibility x °f a strongly correlated Fermi system that 
experiences a rearrangement of the Landau state at crit- 
ical temperature T/. It is instructive to begin on the 
"metallic" side of the transition, where standard Landau 
theory still holds. Within the quasiparticle picture^ 



X 



Xo 



1 - g U 



(4) 



where go is the zeroth harmonic of the Landau spin-spin 
interaction and Xo = ~ f^B^-o is the s P m susceptibility 
of the noninteracting quasiparticles in terms of the Bohr 
magneton ps- The corresponding polarization operator 
is 



o( ) = J -^pT dT = VFPLl{T) 



(5) 



where v F — de^/dp is the quasiparticle group ve- 
locity and Pl is the conventional T- independent, D- 
dimensional Landau result for IIo(T), while 



7(T) = -/ 



dn(Z,T)dp(Q 



(6) 



The derivation of Eq. (0J takes account of the fact that 
the renormalization factor z specifying the quasiparticle 
weight enters only in the phenomenological parameter 
go. According to Migdal's relation T(er; k — 0, ui — > 0) — 



dG- 1 (e)/de 



stemming from spin conservation, ef- 



fects of renomalization of the vertex parts T and the 
Green functions G cancel. 

Noting that dn(£,T)/ 'd£ decays rapidly at £ > T, 
Rolle's theorem may be invoked to reduce integral © 
to the value of the function dp(£)/d£ at a certain point 
£n ~ T. In strongly correlated Fermi systems that 
exhibit non-Fermi-liquid behavior, there are two dis- 
tinct temperature-density regions in which the func- 
tion (dp(£) / d£) £ R shows qualitatively different departures 
from the Landau norm. The first region is situated close 
to the critical density poo. A portion of spectrum 
adjacent to the Fermi surface is flattened anomalously, 
and the curve dp(£)/d£ bends downward, producing In- 
dependent corrections to I(T), and hence to x(^); that 
become important at sufficiently high T. In 2D liquid 
3 He, this behavior is successfully described by the phe- 
nomenological formula21 X (T) ~ [(T**(p)) 2 + T 2 ] ~ 1/2 . 
Still within the first region, Landau theory fails at low 
T ~ as well. New terms in the Taylor expansion, 



Pf , & (p-Pf) 



M*{p) 



■i 



P F 



+ ■ 



(7) 



odd in p — pf, come into play since the first term is 
suppressed by the divergence of the effective mass. Ne- 
glecting the M* term, the cubic term becomes domi- 
nant and upon inserting Eq. (JJJ into Eq. © one finds 




FIG. 1: Derivative dp(£)/d£, in units of (u^) -1 , plotted as a 
function of £/s%, based on calculations in the Nozieres model 
at different temperatures T. Open circles label the points (£r, 
dp/d£\n) on the corresponding curves. Inset: The integral 
J(T), evaluated as (dp/d^)^ R , is drawn versus T/e° F . 



n (T 0) T -2 / 3 . In the second region, T < T f (p) 

and the FC has come to occupy a significant fraction of 
the Fermi sphere. In this case, substituting Eq. © into 
Eq. © yields n (T -► 0) l/T. 

The foregoing inferences are confirmed in numerical 
calculations. For example, Fig. ^ shows results for the 
functions dp/dt; and I(T) based on the Nozieres model^ 
in which the effective interaction between particles is pro- 
portional to S(pi — P2). The transitional behavior near 
Poo is absent in this model, and the function I(T) is sim- 
ply the sum of a term ~ l/T and a constant term. 

In arriving at these results, we have neglected damp- 
ing, which is more pronounced in systems containing a 
FC than in ordinary Fermi liquids. However, in the en- 
ergy region e ~ T relevant to the current problem, the 
damping parameter is moderate, with j{T) ~ T. As 
a consequence, damping does not affect the predicted 
form of no(T — > 0), instead altering only numerical 
factorsii^iS In the transitional region where is not 
completely flat, the ratio -f(T)/T is below 1; hence the 
inclusion of damping has little effect. 

We conclude that there exist two regimes in which 
the spin susceptibility shows characteristic variation with 
temperature, in addition to the standard Fermi-liquid 
domain. The first is a transitional region close to the 
point of fermion condensation, where the FC is incip- 
ient or still minute. Upon inserting the above result 

n (T 0) T~ 2 / 3 into the Eq. Q, simple algebra 

leads to 



X (T -> 0) ~ (C t T 2 / 3 - e t )- 



(8) 



where 0t ~ —gop. Such a power-law T-dependence is 
observed in a number of heavy-fermion metals^*2» 

As the density increases, we pass from the transitional 
regime characterized by iJSJ to the region where the FC 
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FIG. 2: Magnetization for 3 He fi lms at different densities. 
Experimental data from Refs. 1 2 are indicated by symbols, 
solid curves show the present predictions at low T, and dashed 
lines indicate the phenomenological fit of Ref. |21J at higher T. 



contribution plays the dominant role. On combining 

Eqs. 0, ©, and the result U a (T ->• 0) l/T, we 

obtain the Curie- Weiss formula 

X (T -> 0) ~ (T - e^)" 1 , (9) 

with a Weiss temperature ~ ~9o(p)p- Remarkably, 
working within the general framework of the Landau 
quasiparticle approach, we have found conditions under 
which a strongly correlated Fermi liquid behaves as a gas 
of localized spins. 

As a rule, the value of the Weiss temperature @w is 
small (e.g., Qyy — — 0.3mK for the uppermost curve 
in Fig. [21). The reason for this is that there is no 
ferromagnetic instability in any of the systems under 
consideration, i.e. the Pomeranchuck stability condition 
l+goN(0) > is not violated. Since the density of states 
N(0) at the Fermi surface is proportional to M* , preser- 
vation of this condition implies that go(p) changes its 
sign at the critical density Poo where M* — > oo. Hence, 
at this density the character of the correlations changes 
from ferromagnetic to antiferromagnetic; such a change 
does take place in strongly correlated 2D liquid 3 He (see 
Fig. 1 of Ref. Q). 

The predicted behaviors © and © explain the mea- 
surements that trace the gradual evolution of the spin 
susceptibility of a 2D electron layer in silicon as the den- 
sity p is lowered^ The constant susceptibility of the ordi- 
nary Fermi liquid gives way to the T-dependent behavior 
x(T) ~ T _1 of a gas of localized spins. The same trends 
are seen in the extensive data on the magnetization of 
3 He films reported in Refs. [112, some of which are plot- 
ted in Fig. |21 along with our theoretical predictions. 

Before focusing on these results, it must be remarked 
that comparison of different sets of data on 2D liquid 
3 He reveals a crucial dependence of the effective-mass 
enhancement on the number of 4 He layers separating 
the 3 He monolayer from the graphite substrate 2 ^: The 
smaller the number, the greater is the increase of M* . It 



would appear that the enhancement of M* is somehow 
associated with the crystal lattice of the substrate. A 
similar situation occurs in the flattening of single-particle 
spectra of 2D superconductors. For example, in the case 
of a quadratic lattice, flattening is most evident in the 
region adjacent to the van Hove pointsiii*i£ As seen in 
ARPES data 2 ^ on Na . 5 VO 2 , flat portions of the single- 
particle spectrum exist in 2D metals having the same 
Brillouin zone as 2D liquid 3 He. Therefore we suggest 
that the formation of FC in certain segments of the Fermi 
line of 2D liquid 3 He and, correspondingly, the emergence 
of Curie- Weiss components in its spin susceptibility, be- 
gins at densities markedly lower than the threshold for 
fermion condensation in homogeneous 2D liquid 3 He. 

Turning now to Fig. [2 we first consider the uppermost 
curve, belonging to the highest density p = 0.058 A -2 , 
where it is assumed that the FC occupies a significant 
part of the Fermi sphere. The corresponding contribu- 
tion to the spin susceptibility %(T — > 0) is evaluated by 
means of Eq. @. To accomodate higher T ~ Tt, we sup- 
pose that upon approach to Tf the FC density behaves 
as (1 - T 2 /?/) 1 / 2 , with T f = 30mK. (In fact, the spe- 
cific form of this attenuation factor is immaterial.) The 
contribution to TLq(T) from the remaining, normal part 
of the Fermi liquid is assumed to be constant at a value 
~ 0.92. Given these two parameter choices, the experi- 
mental data for x(T,p — 0.058 A -2 ) are well reproduced 
over more than two orders of magnitude on the mK scale. 

The other two curves in Fig. |2 for the densities p = 
0.052A- 2 and 0.036A" 2 , are determined under the as- 
sumption that FC is poised on the verge of onset, so that 
the cubic term prevails in the series representation (J7J| of 
the spectrum £(p). As for the uppermost curve, the effect 
of flattening on x(T — > 0) in this transition regime is at- 
tenuated with rising T through a factor (1 — T 2 /! 12 ) 1 / 2 . 
The constant contribution to Ho(T) from other effects 
is taken as ~ 0.04 for p = 0.052 A" 2 , and ~ 0.014 
for p = 0.036 A -2 . The remaining parameter choices 
are: T t = 20 mK, C t = 1.2 mK 1 / 3 for p = 0.052 A" 2 ; 
T t = 10 mK, C t = 17 mK 1 / 3 for p = 0.036 A" 2 ; and 
t — at both densities. 

When the temperature is lowered to values such that 
T < Tz = pbB, where B is the strength of a magnetic 
field applied to the sample, the perturbative approach 
employed above in evaluating xCO becomes question- 
able. One then needs to solve equations for the spectra 
£±(p) directly. These equations may be written in the 
form 

£+(p) = $ + \t z + J /(p )Pl )i[n+(pi)+n_(pi)]dri , 

£-(p) = Z+{P)-Tz , (10) 

where n + and n_ are the quasiparticle momentum dis- 
tributions of the subsystems with spin projections 
The interaction function / is assumed to be independent 
of the momentum distributions n + and n_ themselves. 
If the field B is small, the solution of Eqs. (|10fl can be 
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FIG. 3: Single-particle spectra £+(p) (dashed line) and £-(p) 
(solid line) in units of £%, together with momentum distribu- 
tions n+ (p) (dashed line) and n_(p) (solid line), plotted as 
functions of p/p F for T z /e° F = 10~ 2 and T/e° F = 10" 3 . The 
Nozieres model is assumed. 



expressed in terms of the solution of the FC problem 

d(p)=$ + J /(p.Pi)no(Pi)dri (11) 

in the absence of the held. Setting B = in Eqs. 1)10(1 . we 
are led to the conclusion that the sum [n+(p) + ri-(p)) /2 
must coincide with the solution no(p) of Eq. f| 1 1 p . 
But according to the second of Eqs. (fTTff) . only one of 
the two states with the same momentum p can be- 
long to the FC. In particular, if no(p) > 1/2, then 
n-(p) = 1 and n + {p) — 2no(p) — 1; hence, in the 
FC region one obtains £-(p) — — T z , while — 
Tln[(l-no(p))/(no(p)-l/2)]. Conversely, if n Q (p) < 
1/2, then n-(p) = 2n (p) and n + (p) — 0; hence, in 



the FC region one has £_(p) ~ Tln[(l/2-n (p))/n (p)], 
while C+(p) — Tz- Thus we infer that at the point where 
no(p) = 1/2, there is a jump in the single-particle ener- 
gies. This behavior is illustrated in Fig. [3 where the par- 
ticle energies £+ (p) and £_ (p) are plotted for the Nozieres 
model with coupling constant f/e% = 0.2. 

Having obtained these results, the FC magnetization 
A4 is readily evaluated as sum of the integral over Uq{p) 
from no(p) = to 1/2, and the integral of 1 — no(p) 
from tiq(p) = 1/2 to 1. Thus we see that in the case 
fifiB > Tz, the FC magnetization M* is independent of 
both T and B. This property implies that the spin sus- 
ceptibility, determined as the derivative \ = dAi/dB, 
does not contain the FC contribution. We therefore con- 
clude that at finite values of B and low T, the spin sus- 
ceptibility x(T) of a system containing a FC increases 
with temperature until T attains values comparable to 
Tz, and then begins to fall off as T _1 if the FC fraction 
is significant. 

In summary, we have shown that if the interaction be- 
tween particles is strong enough to induce a rearrange- 
ment of the Landau state, then the behavior of the spin 
susceptibility \(T — + 0) progressively changes from that 
inherent in Landau Fermi-liquid theory to the behavior 
characteristic of a gas of localized spins. This finding 
calls for reexamination of many aspects of the theory of 
metallic oxides and heavy metals with localized spins. 
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